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There is a simple and general experimental protocol to generate slow granular flows that exhibit 
wide shear zones, qualitatively different from the narrow shear bands that are usually observed in 
granular materials . The essence is to drive the granular medium not from the sidewalls, but to split 
the bottom of the container that supports the grains in two parts and slide these parts past each 
other. Here we review the main features of granular flows in such split-bottom geometries. 

PACS numbers: 45.70.Mg, 47.57.Gc, 83.50.Ax, 



Granular media exhibit a complex mixture of solid 
and fluid-like behavior, often hard to predict or capture 
in models. Perhaps the most striking feature of gran- 
ular flows is their tendency to localize in narrow shear 
bands [1 . A decent model of grain flows should be able 
to capture, and preferably, predict this type of behavior, 
but at present there is no general approach which, for 
given geometry, driving strength and grain properties, 
predicts the ensuing flow fields. 

In recent years, much progress has been made for fast 
flows, such as avalanche flows down an incline, where 
large flowing zones form. Microscopically, momentum 
exchange then takes place by a mixture of collisions and 
enduring contacts. This allows the definition of a dimen- 
sionless parameter /, the inertial number, which charac- 
terizes the local "rapidity" of the flow. A local relation 
between stresses, strain-rates and / then successfully cap- 
tures many aspects of these rapid granular flows [2H4] . 

In contrast, the situation for slow flows, such as those 
made by slowly shearing the boundaries of a container 
containing grains, is still wide open. The averaged 
stresses and flow profiles become essentially independent 
of the flow rate, so that constitutive relations based on 
relating stresses and strain rates are unlikely to capture 
the full physics. In this regime, shear banding is generally 
very strong, with shear bands having a typical thickness 
of five to ten grain diameters. These shear bands often 
localize near the moving boundary. For a recent review, 
see pQ. Experimental handles for probing this shear lo- 
calization appears to be limited, since shear banding ap- 
pears so robust. For example, granular flows in Couette 
cells always show the formation of a narrow shear band 
near the inner cylinder, irrespective of dimensionality, 
driving rate, or details of the geometry [5]. 

In this regime of slow flows, the inertial number / tends 
to zero and momentum transfer is dominated by enduring 
contacts. Soil mechanics is then a natural starting point 
to describe these flows, and both rate independence and 
shear banding are consistent with a Mohr- Coulomb pic- 
ture where the friction laws acting at the grain scale are 
translated to the stresses acting at coarse-grained level. 
The idea is that when the ratio of shear to normal stresses 
is below the yielding threshold, grains remain quiescent, 
while in slowly flowing regions the shear stresses will be 
given by a (lower) dynamical yield stress. This way of 



thinking readily captures the maximal slope of dry sand 
piles. However, the steep gradients associated with nar- 
row shear bands are difficult to capture by a continuum 
theory, and shear bands often are described as having 
zero width [6]. 

Shearbands, then, are not always narrow. In this pa- 
per we will review recent experiments, numerical work 
and theoretical descriptions of wide shear zones which 
have been generated in so-called split-bottom geometries. 
The essence is to drive the granular medium not from 
the sidewalls, but to split the bottom of the container 
that supports the grains in two parts that slide past each 
other. By taking advantage of gravity to drive the granu- 
late from the sliding discontinuity in the bottom support 
of the grain layer, one effectively pins a wide shear zone 
away from the sidewalls. The resulting grain flows are 
smooth and robust, with both velocity profiles and the 
location of the shear zones exhibiting simple, grain inde- 
pendent properties. 

The outline of this paper is as follows. In section [T] we 
will review the results of recent experiments and numeri- 
cal work on the flows which have been generated in these 
special flow geometries. In section [II] we will discuss the 
theoretical ideas proposed to capture the flows observed 
in these geometries. 



I. SLOW FLOWS IN THE SPLIT-BOTTOM 
GEOMETRY: PHENOMENOLOGY 

A. General Description 

In this section, we focus on the rate independent 
regime which is reached for slow enough driving. Two 
variants of the split-bottom geometry will be encoun- 
tered: in experiments one typically employs a cylindri- 
cal split-bottom shear cell, consisting of a bucket, at the 
bottom of which a disc rotates with respect to the bucket 
(Fig. [T^t) [TUTU] , while for theoretical studies the linear 
split-bottom cell with periodic boundary conditions is 
more convenient (Fig. [TJd) [TTHT4] . 
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FIG. 1: (a) Cylindrical split-bottom geometry, showing a disc 
of radius R s at the bottom of a granular layer of depth H. 
Here, the outer cylinder rotates with rate Q and the bot- 
tom disc is kept fixed - A similar geometry, with fixed outer 
cylinder and rotating disc is also frequently encountered, (b) 
Linear split-bottom geometry, where a container is split along 
a straight line in its bottom. This geometry can be seen as 
the R s — » oo limit of the cylindrical cell, (c) The transition 
in flow structure from shallow to deep flows in the cylindrical 
split-bottom geometry. In the dark grey region the material 
essentially co-moves with the disk. 



B. Parameters and Regimes 

The cylindrical split-bottom geometry is characterized 
by three parameters. The radius of the bottom disc R s 
and its rotation rate Q are generally fixed in a set of 
experiments, and the relative motion of disc with respect 
to the cylindrical container drives the flow. The thickness 
of the granular layer, H, is the control parameter that 
typically is scanned in a series of experiments. Note that 
the radius of the container appears immaterial, as long 
as it is sufficiently large; 25% larger than R s appears to 
be sufficient [5]. 

We denote the ratio of the averaged azimuthal velocity 
of the grains vo(r)/r and the disk rotation speed fl by 
u; uj = thus corresponding to stationary grains, while 
u = 1 corresponds to grains co-moving with the driving. 
For the small Q of interest here (typically less than 0.1 
s _1 ), the flow profiles cj(r, z) are independent of Q - the 
flow is rate independent, and transients are short lived. 
Since centrifugal forces are negligible for typical sizes of 
R s (typically a few cm), one may also fix the disc and 
rotate the bucket, and essentially obtain the same sort 
of flows with £j(r,z) = 1 — cj(r, z) - see [THS]. The two 
parameters H and R s set the large-scale structure of the 
flow. 

When the disk rotates, a shear zone propagates from 
the slip position R s upwards and inwards. The qual- 
itative flow behavior is governed by the ratio H/R S1 
and three regimes can be distinguished - Fig. A 
regime of shallow layers is found for H/R s < 0.45, and 




r (mm) 



FIG. 2: (a) Surface flows for glass beads of diameter 300 /xm 
and a range of filling heights H as indicated are well described 
by an error function (fit not shown) - R s = 85 mm here, and 
the outer cylinder is rotating, (b) Collapse of the surface flow 
profiles shown in (a) and comparison to error function. The 
rescaled radial coordinate A is defined as (r — R c )/W. Top 
inset: strain rates are Gaussian. Bottom inset: the tail of 
flow profile corresponds well to the Gaussian tail of the error 
function. Figure adapted from Ref. [7]. 



here the shear zone reaches the free surface. The three- 
dimensional shape of the shear zones is roughly that of 
the cone of a trumpet, with the front of the trumpet 
buried upside down in the sand. Another regime of deep 
layers plays a role for H/R s > 0.65, and here the shear 
zone essentially forms a dome-like structure in the bulk of 
the material; little or no shear is observed at the free sur- 
face. In between there is an intermediate regime, where 
the shear in the bulk of the material is a mix between 
the trumpet and dome-like shape. 



C. Surface Flow 

Shallow layers — We first focus on the flow ob- 
served at the free surface. For shallow layers, a narrow 
shear zone develops above the split at R s , and when H is 
increased, the shear zone observed at the surface broad- 
ens continuously and without any apparent bound. Ad- 
ditionally, with increasing H, the shear zone shifts away 
from R s towards the center of the shear cell (Fig. 2a). 

After proper rescaling, all bulk profiles collapse on a 
universal curve which is extremely well fitted by an error 
function: 
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where erf denotes the error function, r is the radial 
coordinate, R c the center of the shearband (where uj(r) = 
0.5) and W the width of the shearband (Fig.|2|. Accurate 
measurement of the tail of the velocity profile further 
validate Eq. 0, and rule out an exponential tail of the 
velocity profile (Fig. [2|d). The strain rate is therefore 
Gaussian, and the shear zones are completely determined 
by their centers R c and widths W. 

Particle shape does not influence the functional form of 
the velocity profiles [8] , in contrast to the particle depen- 
dence found for wall-localized shear bands in a Couette 
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FIG. 3: (a) Shear zone positions R c versus if, where R s = 
95, 85, 75, 65, and 45 mm. Lines correspond to Eq. [5] (b) 
Log-log plots of W for spherical glass beads of increasing sizes 
(ranging from average diameter 300 /im to 2 mm) for R s = 95 
mm. The lines shown in (b) correspond to exponents of 1/2, 
2/3, 1. Figure adapted from [5]. 



cell [5]. For these, the vicinity of the wall induces par- 
ticle layering, in particular for monodisperse mixtures. 
Apparently such layering effects play no role for these 
bulk shear zones, where it should be noted that the ef- 
fects of shear-induced ordering of particles with larger 
aspect ratios has not yet been investigated. 

Remarkably, the center of the shear zone, R C1 turns out 
to be independent of the material used (Fig. [3^i). There- 
fore, the only relevant length-scales for R c appear to be 
H and R s . We find that the dimensionless 'displacement' 
of the shear zone, (R s — R c )/R s , is a function of the di- 
mensionless height (H/R s ) only. The simple relation 



(R s — R c )/R s — (H/R, 



N5/2 



(2) 



fits the data well (Fig. [3p 

The relevant length scale for the shear zone width W 
defined above is given by the grain properties, and is 
independent of R s (Fig.^). Grains shape, size, and type 
also influence W(H): irregular particles display smaller 
shear zones than spherical ones of similar diameter. The 
best available experimental data shows that 



W/d ~ (H/d) 



2/3 



(3) 



where d denotes the particle diameter. Although this 
scaling has not been checked over more than a decade, 
the exponent is clearly not equal to 1/2 or one. As yet 
there is no explanation for this scaling. 

The available numerical data coming from molecular 
dynamics simulations essentially confirm this picture 
[ID] \12\ \13\ [T5] : the surface flows in split-bottom 
geometries for H/R s < 0.45 are given by Eqs. [l][3j Only 
the absolute width of the shear zone at the surface 
W(z = H) remains as a fit parameter, but once this 
width has been measured for a single value of H/R S1 
Eq. [3] can be used to estimate the width for the whole 
range of H/R s < 0.45. 



FIG. 4: Core precession in a split-bottomed geometry, (a-d) 
Series of snapshots of top views of a setup with stationary 
disc and rotating outer cylinder (for R s = 95 mm, H = 60 
mm, and rotation rate Q = 0.024 rps), where colored particles 
sprinkled on the surface illustrate the core precession for t = 
s (a), t = 10 s (b), t = 100 s (c) and t = 1000 s (d). (e) Data 
collapse of the precession rate uj p for R s — 45 mm (diamonds), 
R s — 65 mm (x) and R s — 95 mm (circles) when plotted as 
a function of H/R s . Figure adapted from Ref. [9]. 



Deep Layers — When H/R s is small, the core ma- 
terial rests on and co-moves with the center disc. With 
increasing H/R S1 the width of the shear zone grows con- 
tinuously, and its location moves inward towards the cen- 
tral region (Eq. [2J. This implies that for deep layers 
qualitatively different flow patterns can be expected to 
occur. 

The most striking feature of these flows is that the core, 
as observed at the free surface, precesses with respect to 
the bottom disc for H/R s > 0.65 — hence material in 
the central part of the system no longer rests on the disc, 
and there is torsional failure of the core. Precession is not 
simply a consequence of the overlap of two opposing shear 
zones, since before being eroded by shear, the inner core 
rotates as a solid blob for an appreciable time (Fig.|3^-d). 

The precession rate u p is defined as the limit of u(r) 
for r going to zero, where we assume, for simplicity, that 
the outer bucket rotates with rate £1 and the bottom 
disc is kept fixed, as in [9] - consistent results are found 
in a setup where the disc was rotated and the outer 
cylinder kept fixed [ID] . For various slip radii, the onset 
height for precession grows with R s , and the data for 
u p collapses when plotted as a function of H/Rs, see 
Fig. [4^). When H/R s becomes of order one, the whole 
surface rotates rigidly with the rotating drum and all 
shear takes place in the bulk of the material; on the 
other hand, for H/R s < 0.65, hardly any precession can 
be observed. 



Intermediate regime — In the intermediate regime, 
0.45 < H/R s < 0.65, a precursor to the transition to 
precession can be observed, since a careful analysis re- 
veals that the surface velocity profiles uo(r) increasingly 
become asymmetric for H/r s > 0.45. In fact, one can 
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FIG. 5: Surface velocity profiles uj(t) for R s = 95 mm and 
increasing layer depth H. Thick curves: H = 10, 20 ... 80 
mm; Thin curves H — 15, 25, ... 75 mm; Dashed curves H — 
56, 57, . . . 69 mm. (a) Precession gradually sets in for H > 60 
mm. (b) Corresponding profiles of x( r ) (dots, see Eq. 
compared to cubic fits given by Eq. [5] (curves). Figure adapted 
from 0. 



generalize Eq. [T] by writing 



^0) = o + ^r/(x(r)) 



(4) 



and by fitting the data for uj(r) over the whole range of 
H/R s to this equation (Fig. 5), one finds that x( r ) can 
be fitted well by a cubic polynomial of the form 



X(r) = a + a 1 r + a 3 r 3 



(5) 



For shallow layers, as = 0, and ao and a\ follow from 
the scaling laws Eqs. [2]and[3j For 0.45 < H/R s < 0.65, 
as starts to grow and governs the symmetry breaking of 
the flow profiles, while for deep layers (H/R s > 0.65), a\ 
tends to zero, and a two parameter fit describes the flow 
profiles well again [9]. 



D. Bulk Flow 

Shallow flows — The bulk structure of granular flows 
is harder to access, but by now, we have information on 
split-bottom flows from experiments that bury and ex- 
cavate colored beads [8 , Magnetic Resonance Imaging 
(MRI) [101 HB] and numerical simulations [TO ) [T2 | [T3 | [T5] 
(see Fig. 6). First, for shallow layers, the flow profiles 
at fixed depth z below the surface H still takes an er- 
ror function form, which allows us to characterize uj(r) 
at fixed z with the same two parameters R c and W as 
before. 

The location of the shear zones in the bulk where found 
to be consistent with a scaling argument put forward by 
Unger et al. The idea is follows: Eq. [2] gives the location 
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FIG. 6: Contours of constant angular velocity uj, for different 
filling height H. Upper panels: MRI experiments: uj — 0.84 
(diamonds), 0.24 (squares), 2.4 x 10~ 2 (circles), 2.4 x 10~ 3 
(triangles), and 2.4 xlO -4 (inversed triangles). Dashed lines 
indicate H and dotted lines are guides to the eye. Lower pan- 
els: simulations. Color is used to identify velocity ranges: 
dark red, uj G [0.84,1]; orange, uj G [0.24,0.84]; yellow, 
uj G [2.4 x 10" 2 ,0.24]; green, uj G [2.4 x 10" 3 ,2.4 x 10" 2 ]; 
blue, uj G [2.4 x 10" 4 ,2.4 x 10" 3 ]; grey, uj G [0,2.4 x 10" 4 ]. 
Figure reprinted with permission from [10]. Copyright (2006) 
American Physical Society. 



R c at the free surface. Then to obtain R c at depth z, 
one imagines a systems with a depth of H — z, and by 
inverting Eq. [2j obtains where the split would have to be 
in a system of depth of H — z for the surface location to 
be as given [6] . Identifying the split at depth H — z with 
the center of the shearband R S (H — z), this yields: 



H - R c [1 - R c /R s (l - H/R s ) 2 - 5 )] 



1/2.5 



(6) 



The width W(z) of the shear zones is harder to obtain 
reliably, but the best available experimental data suggest 
a power law growth of the form W ~ z a , where a is 
less than 1/2 and more than 1/4 [SJ [11]. More recent 
numerical studies ^3 found that W(z) can also be well 
described by a "quarter circle" curve of the form 



W(z) = W(z = H)^l - (1-z/H) 2 . 



(7) 



Deep Layers — The symmetry breaking and the even- 
tual disappearance of grain motion at the surface indi- 
cates that qualitatively different bulk flow is developing: 
the trumpet shape of the shear zones present in shallow 
layers must have changed. When H/R s is sufficiently 
large, the shear zone is entirely confined to the bulk of 
the material, and forms a dome-like structure above the 
rotating disk [6j EJ [TO] US] - see Fig. [6] The torsional 
failure of the material is thus concentrated in the dome. 
Cheng et. al. measured the functional form of the axial 
velocity profile uj(z)\ r= o : and found it to be described by 
a Gaussian: 



uj(z, r — 0) = uj p + (1 — ojp) exp — z 2 / (2a 2 ) 



(8) 



where uj p is the rotation rate observed at the surface at 
r = 0, which decreases roughly exponentially with i7, 
and <j, the width of the bulk Gaussian velocity profile, is 
approximately R s /h [10]. More recent experiments find 
a slightly different functional dependence on z [T7] . 

Couette versus split-bottom geometries — 
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FIG. 7: Surface flow profiles observed for 1 mm glass beads 
in a split-bottom Couette geometry, with inner cylinder ra- 
dius Ri — 40 mm, a split at R s — 65 mm, and H — 
10,30,40,50,60,70,80,100,110 mm. The outer cylinder is 
120 mm. Figure from Ref. [18]. 



The first studies [7] of split-bottom geometries where 
done in a modified Couette cell, as shown in Fig. [7| The 
resulting flows are similar to the disc geometry, as long 
as the shear flow is away from the side walls, but since 
for increasing filling heights the shear zones move inward, 
they will inevitably 'collide' with the inner cylinder for 
sufficiently large filling height. The resulting flow profiles 
are shown in Fig. [7] [18] • 

First, one observes that for sufficiently large H, the 
flow profiles become independent of H. The main result 
is that the tail of these flow profiles becomes purely 
exponential for large H, while it is Gaussian for shallow 
H. We have found this exponential tail to be robust, 
i.e., independent of grain shape. Note that this does not 
contradict the findings of Mueth et al - these concern 
the shape of the flow profile near the shearing wall, 
corresponding to the range 10 -3 < uo < 1 - flow profiles 
are indeed grain dependent here. But further out in the 
tail they become purely exponential. For other examples 
of exponential tails see [19] . 
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FIG. 8: Evolution of dilatancy. (a) Color map of relative 
density change (light blue corresponds to -10%) for H/R s = 
0.51 after 4 rotations of the bottom disk (b) Spread of dilated 
zone for vertical shear observed in the bulk at H/R s — 0.51 at 
z = 22 mm (H 36 mm), for N = 1/2, 1, 2, ... , 64, compared 
to estimates where, after N turns, the local strain equals one. 
(c) Spread of dilatancy for dome-like flow observed at H/R s = 
0.77, for number of disc rotations N as indicated. Figure 
adapted from Ref. [TB] . 



the strain rate [20] . 

Consistent with this, the dilated zone was found to 
slowly spread throughout the system as time progresses 
(Fig. [8]). This spread is consistent with the idea that, af- 
ter initial preparation, the accumulated local strain gov- 
erns the amount of dilatancy. To show this, the flow field 
in the bulk was reconstructed by combining the above 
mentioned scaling relations for the location and width of 
the shear zones in the bulk, and this reconstructed flow 
field can then be compared to the density field obtained 
by MRI. Such comparison shows that the locations of the 
dilated zone and the shear zone coincide, for small filling 
heights (H/R s < 0.6). Finally, for deep filling heights 
where torsional failure and precession play a role, a rela- 
tively long-lived transient was found to cause the dilated 
zone to deviate substantially from the late-time shear 
zone. 



E. Dilatancy 

By means of MRI, direct measurements of the evolu- 
tion of the local packing density of the shear flow gener- 
ated in a cylindrical split-bottomed geometry were per- 
formed in [16 . To be able to image the particles by means 
of MRI, food grade poppy seeds were used; these contain 
MRI-detectable mineral oils. 

It was observed that the relative change in density in 
the flowing zone is rather strong and saturates around 10- 
15 % - likely the rough and peanut shape of the poppy 
seeds plays a role. After long times, a large zone with 
essentially constant low packing fraction forms, encom- 
passing most of the shearband. The fact that the density 
remains constant here, even though local strain rates vary 
over many decades, suggests that the density of the flow- 
ing material is a function of the total strain, and not of 



II. FLOWS IN THE SPLIT-BOTTOM 
GEOMETRY: THEORY 

The flows in split-bottom geometries have been 
simulated both by molecular dynamics simulations 
[T0l [T2| [T5] as well as by contact dynamics [13 , and 
a number of theoretical approaches have been put 
forward. It remains remarkable that no single consistent 
theoretical framework is available from which Eqs. [T]3] 
can be deduced. Creating such a theory would constitute 
an important milestone in the development of our un- 
derstanding of slow granular flows. Here we discuss the 
main approaches to flows in the split-bottom geometry. 

Variational principle — The first attempt to de- 
scribe the flows in split-bottomed geometries goes back 
to Unger and coworkers [6]. The flows are treated in a 
Mohr-Coulomb fashion, with shear bands of zero width. 
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The idea is to minimize the energy dissipation needed 
to sustain the flow. Calculating the total friction along 
the shear-sheet r(z) using assumptions of constant fric- 
tion coefficient \i and hydrostatic pressure P, amounts to 
finding the minimum of the functional (see Fig. 9) 



T(H) = 2gir P (j)n / (H - z)r 2 \/l + (dr / dz) 2 dz (9) 
Jo 

Here p is the bulk density of the particles, <p is the av- 
erage packing fraction (~ 0.59 [4 j) and \i is the effective 
friction coefficient. Identifying r(z) with the center of 
the shearband R c (z), minimizing Eq. [9] gives predictions 
for the location of the shearbands in the split-bottom 
geometry. The location of the shear zones predicted for 
shallow layers are very close to those observed, and for 
deep layers the model predicts a hysteretic transition to 
dome-like shear. Hence, while a number of aspects of 
split-bottom flows can be captured by this simple model, 
the hysteretic transition and zero width of the shear 
bands are clearly in contrast to experimental findings. 

Extending the variational principle — To ex- 
plain the broad shearbands in the split-bottom geome- 
try, a random or randomly varying local material failure 
strength p~4j [21] was invoked. The main extension ex- 
tends the minimal dissipation model by combining the 
variational principle with a self organized random po- 
tential as follows. At any given time, the shear band is 
represented as having zero width. However, the granu- 
lar material is now taken to be inhomogeneous, with a 
local strength field which varies with space and is up- 
dated every timestep. At any given time the surface 
that minimizes the torque can be obtained, after which 
the strength field is updated etc. A smooth flow pro- 
file is then obtained by averaging over the different shear 
bands. The resulting flow fields are very close to those 
observed experimentally, with the only adjustable param- 
eter controlling the effect of disorder. One possible point 
of criticism is that the model assumes that the fluctu- 
ating shear bands are radially symmetric - see [14] for 
details. 

The main findings of [21] where confirmed by a 
different but related approach by Jagla [14 , and recently 
a two dimensional model using stochastic interparticle 
forces was shown to be able to also generate shear bands 
of finite width [22] . 

Inertial flows — Jop performed simulations of 
flows in a cylindrical split-bottom setup [23], using 
the inertial number theory, which should be valid for 
faster flows [4 . The location of the shear zones in 
the bulk, the smooth transition to precession and the 
dome flows were all recovered. The width of the shear 
zones was found to scale with driving rate as ^°- 38 , 
and indeed for slow flows the shear zones attain zero 
width. The inertial model therefore does not fully 
capture the physics of the slow split-bottom flows, but 
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FIG. 9: (a) An arbitrary shear zone of zero width r(z) sep- 
arates a rotating inner core and a static outer body. In (b): 
The frictional stress a r e = I^P on the shearing surface can be 
integrated to give the total driving torque necessary to rotate 
the disk in red. The integral is given in Eq. [9] 
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FIG. 10: (a) Shear free sheets in a linear shear geometry, (b) 
Stress components acting in material (c) Shear free sheets in 
curved geometry, (d) Stress ratios in section of linear geom- 
etry, showing that the stresses are of the form of Eq. [To] and 
that the friction coefficient \i — — (723/033 is not completely 
constant. Figure adapted from Ref. |12j . 



it does slightly better than lingers original model [6 
in that it captures the smooth transition to precession. 
The only experimental work so far on rate dependent 
flows in the split-bottom is the work by Corwin [24] . 
He however described flows for which the centripetal 
forces dominate the gravitational and hydrostatic forces, 
which are outside the regime of flow rates studies by Jop. 

Shear free sheets — The form of the stress tensor 
in slow granular flows is a matter of debate. From a soil 
mechanics perspective, there is no reason to assume that 
the principle directions of stress and strain rate align, but 
Depken and coworkers have suggested that once there 
is flow in the system, the principle directions do align. 
The idea is that flow, even at a distance, creates suffi- 
cient amount of agitation that any amount of shear stress 
should lead to flow - in other words, once there is flow, 
there is no clear yielding threshold anymore [TTJ [12] [17] . 
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Coaxiality severely restricts the form of the stress ten- 
sor, and for steady grain flows the flow can be decom- 
posed into so-called shear free sheets (SFS), that slide 
past each other - there is no (average) grain motion 
within the sheets (Fig. 10i-c). In this SFS basis, the 
stress tensor takes the form: 



structure of this model is incapable of capturing the 
observed wide shearbands - the model's only lengthscale 
is the spot size, and therefore cannot yield wide shear- 
bands. Moreover, its use of Mohr- Coulomb plasticity 
theory is in conflict with the observed co-axiality p~2j [13] 
of the stress and strain tensors. 



/ P' \ 

<TSFS = P T . (10) 

V r P) 

The form of the stress tensor is reminiscent of that of 
the inertial model [3 , but with the exception that the 
stresses P and P' do not need to be equal (they are not 
in fact - see Fig. [l0| ), and that the friction coefficient 
fji : = T /P does not depend on flow rate, since the flows 
are rate-independent. However, as Depken showed [11], fi 
cannot be constant if the shear zones have finite width, 
and in fact has to attain a local maximum within the 
shear zone - as subsequent numerical simulations indeed 
found [rang. 

Recent contact dynamics simulations of Ries et al. 
in linear split-bottom cells [13] have confirmed that 
the stress and strain tensors align, so that the stresses 
take the form given by Eq. [To] The alignment also 
occurs in the absence of gravity (to carry out these 
simulations, a mirror of the system was added; see Fig. 7 
in [13 ). Perhaps surprisingly, gravity appears not to 
be important for the understanding of split-bottom flows. 

Spot model — To describe rate-independent flows 
in general, Bazant and coworkers put forward the 
'spot '-model, which is based on the assumption that 
slow, dense granular flows are best described with a 
diffusion of low density regions in the material, called 
spots [25 . The two-dimensional model describes flow 
profiles in chute flow and Couette geometries reasonably 
well. However, recently it was shown [26 that the very 



III. OUTLOOK 

The surprising, robust and universal properties of slow 
granular flows in split-bottom geometries have made the 
split-bottom geometry into a versatile testing ground, not 
only for models of slow flows, but also for experimental 
studies of related flows, of mixtures and segregation [27] , 
of non-local flows [17 , of faster grain flows [24, 28 and 
of suspension flows [29] . 

So far, no single convincing continuum theory to de- 
scribe the wide shear zones generated in split-bottom 
shear cells has been put forward — even though the 
experimental results strongly suggest that these flows 
should be amenable to a continuum description. The 
failure of the Mohr- Coulomb approach to describe the 
internal structure of narrow shear bands might not be 
troublesome, but its failure to describe these much wider 
shear zones is cause for concern. 

Apart from the theoretical challenges and experimental 
use of the split-bottom geometry, the broad shear zones 
occurring in this geometry allow for further experiments 
on slow flows, that are more difficult to realize in nar- 
rowly localized shear bands. Open questions for the fu- 
ture include to understand the microscopic organization 
and velocity fluctuations in these flows, to understand 
the role of interstitial fluid and grain shape, and to ex- 
plore the range of much smaller, but in particular much 
faster driving rates. 
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